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Abstract 



Hadronic tt+tt" atom is studied in the relativistic perturbative approach based 
on the Bethe-Salpeter equation. The general expression for the atom hfetime 
is derived. Lowest-order corrections to the relativistic Deser-type formula for 
the atom lifetime are evaluated within the Chiral Perturbation Theory. 
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Typeset using REVT^ 



I. INTRODUCTION 



The pion-pion scattering amplitude at low energies forms one of the basic building blocks 
in the hierarchy of strong interaction processes. It serves as a useful probe for the investiga- 
tion of the effect of Chiral Symmetry breaking since in the chiral limit the pion interactions 
vanish at the threshold. According to common belief, low-energy interactions of pions are 
described within the Chiral Perturbation Theory (ChPT) [|l],^ which exploits the full con- 
tent of global QCD symmetries. The vrvr nir amplitude in the ChPT is obtained as an 
expansion in quark masses and external pion momenta. The predictions of ChPT are sensi- 
tive to the magnitude of the quark condensate < 0|gg|0 >. In the standard scheme with 
a "large" condensate, the S'-wave vrvr scattering lengths are predicted to be Qq = 0.217 and 

— Oq = 0.258 in units of the inverse charged pion mass ||^. The calculations within the 
Generalized ChPT with a small quark condensate yield Oq = 0.27 0. Despite a significant 
difference between these predictions, both results for the scattering length are compatible 
with the experimental value = 0.26 ±0.5 0]. Consequently, a precise measurement of tttt 
scattering lengths will be an excellent test of the ChPT. In particular, an experiment of that 
sort would provide an important information about the behavior of the quark condensate in 
the chiral limit which in its turn is related to the properties of the QCD gluon vacuum. 

The experimental study of the vrvr scattering process is a very difficult task mainly due 
to the absence of the pion target. Indirect information extracted from the available data for, 
e.g., the process vrA^ — mvN 0, produces large error bars when extrapolated to the two-pion 
threshold. The study of Kg^ decay which is preferable for determining of the parameters 
of TTTT interaction near threshold, is complicated because of a very small branching ratio of 
this process 0] (for the review of the recent status of vrvr experiments see, e.g. Refs. 0). In 
view of the above-mentioned experimental uncertainty in the determination of the nir data 
near threshold, the forthcoming high-precision measurement of the tt^tt" atom lifetime by 
the DIRAC collaboration at CERN (project # PS212) attracts much attention since it will 
allow the direct determination of the difference af] — Oq and thus will provide an excellent 
probe for the predictions of ChPT. The possibility of observation of such atoms was argued 
in Ref. [0. The first estimation of the lifetime of an atom formed by 7t~^ and vr" in the 
ground IS" state ti = 2.9t^]^ x 10^^^ s was given in Refs. 0. The expected high-precision 
experimental data from DIRAC experiment call for a refined theoretical treatment of this 
sort of bound systems. 

Since the characteristic average momentum in hadronic atoms are of an order of a few 
MeV, these systems are highly nonrelativistic. With the use of this fact the nonrelativistic 
Deser formula was derived in Refs. P,p!0|. For the particular case of the tt+tt^ atom the 
formula reads 



_i _ IGtt /2 Am^y/2^^Q ^2^2 



<-air |v&„(0)p (1) 



which relates the lifetime r„ of the atom in the n-excited state to the value of the Coulomb 
wave function (w.f.) of the atom at the origin \E'n(0) and the difference of the S'-wave vrvr 
scattering lengths with the total isospin 1 = and 1 = 2. Note that in (|^) one assumes the 
isospin symmetry when expressing the scattering length for the reaction 'k~^tt~ ir^n^ in 
terms of scattering lengths with a definite total isospin, though the factor Am^ in the r.h.s. 
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of Eq. (I]) indicates the necessity of taking into account the isospin-breaking effects in the 
theoretical description of the decay process. Note also that the Deser-type formulae for the 
energy level displacement and lifetime are now widely used for the theoretical analysis of 
the experimental data for other hydrogenlike bound systems, such as pionic hydrogen ||11|| , 
pionic deuterium [Q, etc. 

It is well known that the Deser-type formulae for hadronic atoms in the nonrelativistic 
scattering theory are valid up to the electromagnetic corrections to hadron scattering pro- 
cesses, if the mass difference between the charged and neutral components of the isotopic 
multiplet is assumed to be of a purely electromagnetic origin. Put differently, the clear-cut 
factorization of strong and electromagnetic interactions in the hadronic atom observables, 
which is explicit in Eq. (|l]), is valid up to (small) electromagnetic effects. In Ref. |jl3 



the regular approach was constructed for the evaluation of these corrections. The ideas of 
this approach have been successfully applied to the study of the properties of a vr+vr" atom 
in Refs. |n,T^ where the coupled-channel scattering problem with (vr+vr") and (tt'^tt^) free 
pairs in the asymptotic states has been considered. The nonrelativistic scattering theory has 
been used for the investigation of hadronic atoms also in Refs. [p!6| , pT| . Note, however, that 



very "narrow" and "deep" phenomenological tttt potentials used in Ref. lead to the in- 
stability in the calculated observables of the bound state with respect to a small variation of 
input parameters. A strong enhancement of the potential in the vicinity of the origin within 



the inverse scattering theory approach, used by the authors of Ref. |T^, stems mainly from 
the particular parameterization of the vrvr scattering phase shift in the high energy domain, 
where theoretical calculations of this quantity cannot be performed. 

To summarize, the lowest-order Deser formula (|1|) and its counterpart for the energy- 
level displacement of an atom due to strong interactions are valid irrespective of the con- 
crete choice of the strong interaction potential between hadrons. The magnitude, and even 
the sign of corrections to it, however, strongly depend on this choice. For the theoreti- 
cal analysis of the high-precision experimental data expected from DIRAC experiment, the 
model-independent evaluation of these corrections is needed, based on the underlying (chi- 
ral) Lagrangian dynamics of hadrons, rather than the nonrelativistic potential picture of 
strong interactions. 

The problem of relativistic description of hadronic atoms is much richer in content than 
the same problem in the nonrelativistic scattering theory formulation. Many new effects 
which were absent, or were mimicked in the potential treatment, now arise naturally from the 
beginning (e.g., the "vacuum polarization" and "finite size" corrections which are borrowed 
in the potential picture from field theory). Namely, the problem of evaluation of the atom 
lifetime on the basis of the underlying strong interaction dynamics was addressed to in Refs. 



18| , p^ . Ref. [0 deals with the so-called retardation correction in the pionium lifetime. 



In Ref. the radiative corrections to vrvr scattering lengths have been evaluated, which 



induce the corresponding 0{a) correction in the pionium lifetime. Recently, the paper by 
H. Jallouli and H. Sazdjian appeared p2| which is aimed at the consistent description of 
the properties of tt+tt^ atom on the basis of the 3D bound-state equation obtained in the 
framework of constraint theory, with the underlying strong vrvr interactions described by 
ChPT. The authors have calculated the corrections to the pionium lifetime coming from 
the mass difference m^± — 771^^0, as well as the corrections from second-order perturbation 



theory and electromagnetic radiative corrections. In paper p3| the correction due to vacuum 
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polarization to the pionium lifetime was calculated. Below we shall present the detailed 



comparison of our results with those given in Refs. |l22| , |23 



Our previous papers [0,0 were aimed at the consistent field-theoretical treatment of 
TT+TT^ atom observables on the basis of the Bethe-Salpeter (BS) equation. Namely, in Ref. 
2^ we derived the relativistic analogue of the Deser formula (|^) for the pionium lifetime 



and evaluated the correction to it, coming from the displacement of the bound-state pole 
by strong interactions, referred hereafter to as "strong correction". In Ref. pSI] we pre- 



sented a systematic perturbative approach based on the BS equation for the calculation 
of hadronic atom observables. In this paper we give a closed expression, containing all 
first-order corrections to the pionium lifetime, and evaluate part of them. Namely, apart 
from the " strong correction" , which is reproduced here, we calculate the correction due to 
the exchange of Coulomb photon ladders, corresponding to the "second-order perturbation 
theory" correction from Ref. [^], and the relativistic correction to the bound-state w.f. 

The purpose of the present paper is to give a detailed description of the systematic 
perturbative approach to the hadronic atom characteristics, based on the BS bound-state 
equation. This approach is by construction free of any double-counting problems. Within 
this approach we collect together and calculate or give an estimate of all lowest-order correc- 
tions to the Deser-type relativistic formula for the pionium lifetime. The underlying strong 
vrvr interactions within our approach are described in ChPT. Consequently, the results of 
present calculations of the lifetime of a tt+tt^ atom form a self-consistent basis for the verifi- 
cation of the predictions of ChPT in the DIRAC experiment. Further, in the present paper 
we discuss in detail the links to other approaches used for the description of a vr+vr" atom 
and, in particular, that from Refs. p^ , |23[] . 

The layout of the present paper is follows: In Sect. II, we present a detailed description 
of the perturbative approach to the tt'^tt' atom characteristics. In this section, we give a 
closed expression for the first-order correction to the pionium lifetime. In Sect. Ill, we give 
an evaluation, term by term, of various first-order corrections to the Deser formula (|T|). Sect. 
IV contains our conclusions. 



II. PERTURBATIVE BETHE-SALPETER APPROACH TO tt+tt ATOM 

The evaluation of corrections to the Deser-type formulae for hadronic atom observables 
cannot be confined solely to the evaluation of corrections to the pionium w.f. or to scattering 
lengths in Eq. (|I]). One has to develop a consistent perturbative scheme for the calculation 
of atom observables, which in the lowest-order approximation should yield the Deser-type 
relations for these quantities. Our approach is based on the field-theoretical BS equation 
with the kernel constructed from the underlying Lagrangian of ChPT. 

Below we shall briefiy discuss the basic ideas and assumptions of the approach which 
employs the following physical picture. The formation of a 'k~^tt~ atom proceeds mainly due 
to the static Coulomb potential whereas strong interactions are responsible for its decay. 
The atom is described by the exact w.f. which obeys the field theoretical BS equation. 
For our purposes we split this kernel into the Coulomb piece and the remainder, the latter 
considered in our scheme as a perturbation. Then, the exact BS w.f. is related to the 
relativistic Coulomb w.f. in the perturbation theory. The crucial point of our approach 
is that with the use of the above relation the observables of an atom (lifetime and binding 
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energy) in every perturbative order are expressed in terms of the Coulomb w.f. In the leading 
order of a perturbative expansion we reproduce the Deser formulae for atom observables . 
The next-to-leading term in the perturbative expansion produces all lowest-order corrections 
to the Deser formulae. 

Let us now pass to the description of the perturbative BS approach to tt+tt^ atom 
observables. The existence of a quasistable 7t~^tt^ atom (r ~ 10^^^ sec) corresponds to the 
bound state pole in the four-point Green function for the transition vr+vr" — > tt+tt" at a 
complex value of the cm. energy = = — iMT. Here M denotes the "mass" of an 
atom, and F stands for the decay width. Hereafter all formulae are restricted to the c.m.s. 
of an atom. 

The BS w.f. of an atom xbs for — » obeys the exact BS equation (Fig. la) 

G^\P;p)XBsiP;p) = J -0^,VBs{P;p,q)XBs{P;q) (2) 

Here G2{P]p) = D{\P + p)D{\P - p) is the two -pion Green function where D{k) stands 
for the dressed pion propagator. Further, Vbs denotes the BS equation kernel, which is a 
sum of all four-point one-particle irreducible diagrams with amputated external legs. 

It is appropriate to "transfer" the self-energy insertions in the charged pion exter- 
nal legs to the r.h.s of the BS equation. This can be easily achieved if one defines 
xiP;p) = Go{P;p)G^\P;p)xBsiP^p) and ViP;p,q) = VBsiP;p,q)G2iP;q)G^\P;q), 
where Gq{P; p) = i[{^P + p) — m^] x i[{\P — p) — m^] is the free two-particle Green 
function and = m^± denotes the charged pion mass. The diagrammatic expansion of 
the new kernel V is given in Fig. lb. It, in addition to the diagrams included in the "true" 
kernel Vssi contains the self-energy diagrams in outgoing external pion legs, i.e., only half 
of the possible insertions in external legs. Note that this property of the new kernel V is 
crucial for proving the gauge invariance of bound-state characteristics, as well as for demon- 
strating the cancellation of infrared singularities in these characteristics (see below). The 
BS equation for the new w.f. x depicted in Fig. Ic is given by 

G-,\P-p)x{P:p) = j -0^,ViP;p,q)x{P;q) (3) 

The kernel V contains the instantaneous Coulomb part Vc which, in a complete analogy 
with the positronium case, is responsible for the formation of the bound state composed of 
7r+ and vr". We are willing to develop the perturbative expansion of the atom observables 
in the "remainder" of the potential denoted hj V = V — Vc- For this purpose we give 
first a complete solution of the "unperturbed" problem, with the kernel containing only the 
instantaneous Coulomb part. 

We choose the instantaneous Coulomb part of the potential according to the Barbieri- 



Remiddi prescription [^], to be ||2^ 



Vc{p,q) = {w{p)y/' {wm/\ w{p) = {ml+p'f' (4) 

Note that a particular choice of the Barbieri-Remiddi kernel (§) is only the matter of con- 
venience and the final results are not affected by this choice (below we shall demonstrate 
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this property of the perturbative expansion exphcitly). However, choosing as an initial 
kernel, one can take advantage of the fact that the BS equation with a kernel of this sort is 
exactly solvable, with the properly normalized ground-state solution given by [^] 



iJc{p)=iG,{M*;p)A{w{p)fl^ 



Anam 



2 _|_ ^2 



(5) 



where 7 = ^m.„a, tt^q = 7^ and (M*)^ = m^(4 — a^) is the eigenvalue corresponding to the 
unperturbed ground-state solution. The c.m.s. momentum in the free Green function Go in 
Eq. has the components (M*,0). 

The normalization condition for the Coulomb w.f. reads as 



< iJc\N{M*)\iJc >= 1, N{M^-p,q) = {2nY 6^'\p - q) 



d 



2M* dM* 



Go\M*;p) (6) 



and the scalar product in the momentum space is defined by the integral over (i'^g/(27r)^. 
We shall use this shorthand notation hereafter. 

The exact solution for the Green function corresponding to the nonrelativistic Coulomb 
problem, was given by Schwinger |2^. Using this result, one can obtain the solution for the 
4D Coulomb Green function corresponding to the kernel ( 
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Gc{P'- p, q) = (27r)^5W {p - q)Go{P'; p) + Go(P*; p)Tc{E^- p, q)G,{P*- q) 



Where 



Tc{E*;p; q) = IGivrm^a 
Dip; E*;p,q) = ip-qf 



4E* 



[{p- qf 
2 



+ 



vdpp ^ 



D{p; E-;p,q) 



(7) 



(8) 



E' 



v = a 



AE* 



_ P_ 

1/2 



E"" 
E* 



-)(1-P)^ 

IT ' 



1 

m. 

(p.)2 



Ami 



4m, 



The first and second terms in the r.h.s. of Eq. (§) correspond to the exchange of one and 
multiple Coulomb photons, respectively. In the vicinity of the bound-state pole (P*)^ 
{M*y, t/ — s- 1 and the integral in the r.h.s. of Eq. @ diverges as Jq dp/p. Extracting this 
divergent piece which corresponds to the bound-state pole in the Green function, one can 
write 



r ._, MM*;p)MM^;q) 



(9) 



where in the vicinity of the bound-state pole the regular part of the Coulomb Green function 
takes the form 



GniM*;p,q) = (27r)V^)(p - g)Go(M*; 



+ i{w{p)w{q)) 



1/2 



mq)-S{p)S{q) 



8 d 



M* dM* 



Go{M^;p)Go{M*;q) 
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{p- q) 



+ - -,S{p)S{q)mq) (10) 



dp 1 1 

JO p 4 4 

The solution x{P',p) of the exact BS equation (|^) can be expressed via the unperturbed 
solution ipc{M*\p) by the following limiting procedure ||2^j2r 



<X\=C < ^c\Gc\P*)G{P), {P*y {M*Y, P' ^ M' (11) 

where C denotes the normalization constant. Note that this relation is the relativistic 
generalization of the well-known nonrelativistic formula 

< Xl = lim iv < M -T\ Ij , . (12) 

r)^0(+) E — H + 17] 

which connects eigenvectors of the total Hamiltonian H with the unperturbed eigenvectors 
(see, e.g. 



2S 



The result in Eq. (pjj) depends on details of the limiting procedure. This equation 
makes sense if the quantities (P*)^ — (M*)^ and P^ — are assumed to be the infinitesimal 
variables of equal strength. In Refs. [p4| , ^ we have assumed the prescription (P*)^ = 
(M*)^ + A, P^ = Af^ + A, A ^ 0. Note that we can employ this prescription without 
the loss in generality, since the change of the direction in the ((P*)^, P^) plane along which 
this limiting procedure is performed affects only the normalization constant C. Further, the 
validity of Eq. ([TT|) can be trivially checked by extracting the bound-state pole in G(P) and 
using the BS equation for < ipc\- 

Let us now introduce the relativistic generalization of the projector operator onto the 
states orthogonal to the ground-state solution 

g = l - iV(M*)|^c><^c| (13) 



Then, with the use of the Hilbert identity it is easy to demonstrate that Eq. (|TT|) can be 
rewritten as follows 



< Xl = < x\N{M*)\iJc >< i^cl (I + V\P)[Go\P'') -Vc- QV\P)]-^Q 

- < xIAGo '(P^) -Vc- QV\P)]-'Q (14) 

where AGg ^ = Gq^{P) — Gq^{P*) and the limiting procedure is implicit. In the derivation 
of Eq. ([T^) we have used 

< MGc\P*)[G,\P'') -Vc- QV\P)]-' = (15) 

which stems from the fact that the inverse operator in the l.h.s. of this equation does not 
have a pole in this limit. With the limiting prescription chosen above the normalization 
constant equals p4 
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-C-'=< x\N{M*)\^c> 



(16) 



Equation (H) can be solved with respect to < xl, resulting in 

< Xl = < [1 + V'{P)GvQ] [1 + AG^'GvQr' (17) 

where the operator GyQ obeys the equation 

GvQ = Gr{P'')Q + Gr{P*)QV\P)GvQ (18) 

and Gr stands for the regular part of the Coulomb Green function. It can be easily demon- 
strated that Eq. (p]7| ) can be rewritten as 

< Xl = -C-^ < + (AGo^ - V'{P))GrQ]-' (19) 

Substituting this solution into the complete BS equation (H) and using the BS equation for 
the function \ipc >, we arrive at the final relation 

< i;c\[l + (AGo' - V'iP))GRQ]-\AG^' - K'(P))|V'c >= (20) 

which provides the basis for the perturbative expansion of bound-state observables (the 
energy of the atomic level and decay width). Namely, the only unknown quantity in the 
l.h.s. of Eq. (|20|) is the bound-state total four-momentum P, which enters parametrically 
into this expression. Expanding the l.h.s. of Eq. ( [20| ) in the perturbation theory up to a 
given order, one can then determine bound-state observables with a required accuracy. Note 
also that Eq. ( pOD is a complex equation, and in every perturbative order it provides two 
real equations for determining the energy level shift and decay width. 

Equation (^^, however, still contains the BS kernel, and does not contain the vrvr scatter- 
ing amplitudes. Below we carry out the first-order perturbative calculations and demonstrate 
explicitly that only these amplitudes appear in the final result. For this purpose let us note 
first that the quantity GrQ in Eq. (|20|) is given by formulae similar to (p!OD, with $ and R 
replaced by $ and R, respectively, and 

-R(P) q) = -R(P) q) + 5 + higher orders in a (21) 

(this can be demonstrated by straightforward calculations, using Eqs. (||), (1), (p!OD and 
(0)). We isolate the free part in GrQ by writing GrQ = Go{M*) + 6G. 

Let us now turn to the perturbation kernel V'{P). This potential can be decomposed 
into the following parts: 

1. A purely strong part, which is isotopically invariant. This part survives when electro- 
magnetic interactions are "turned off' the theory. 

2. The part which is responsible for the m^± — m^o electromagnetic mass difference 

3. Remaining electromagnetic effects, including the exchanges of virtual photons. 

Parts 1 and 2 are regarded to be more important for the following reasons. The 
first term includes strong interactions which govern the decay of a pionium. The sec- 
ond term makes this decay kinematically allowed. Consequently, it seems to be natu- 
ral to consider them together, denoting the corresponding potential as V12 = Vi + V2. 
The T-matrix corresponding to summation of the potential V12 in all orders is given 
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by TuiP) = Vi2{P) + Vi2{P)Go{P)Ti2{P). The rest of the potential is referred to as 
V3 = V — V12 and is treated perturbatively. 

We would like to emphasize once more that this splitting is rather convention-dependent 
and is dictated by the convenience considerations. In practice it is convenient to include, 
into parts 1 and 2, as much terms as possible. It is obvious, however, that the final results 
do not depend on the prescription chosen for that splitting. 

We perform the perturbative expansion of the basic equation (^) in V3 and 6G up to 
the first nontrivial order. Meanwhile we expand AGg ^ the Taylor series in the variable 
SM = M - M* and substitute 

Expressing everywhere V12 in terms of T12, we finally arrive at the identity 
= ~2iM*5M- < tpc\Ti2\^c > 
+ 6M< V^cKGo ')'G'oTi2|V^c > +^(5M)2 < V^cKl + T,2Go)iG,')" + {G.YGoTuli^c > 

+ < ilJc\{6M{Go'y - Tu)6GTi2\^c > - < tpc\{l + TuGo)Vs{l + GqTuMc > (23) 

where Go = Go{M*) and the prime stands for the differentiation with respect to M*. 

Equation p3| ) contains all first-order corrections to the pionium lifetime. Below we 
examine this relation term by term. Note that we employ the commonly accepted "local" 
approximation, i.e. we assume that the quantity T12 does not depend on relative momenta. 
The origin of this approximation can be traced to the "sharpness" of the Coulomb w.f. of an 
atom which in the momentum space has the characteristic range 7 ~ 1 MeV, much smaller 
than the typical hadronic scales. 

III. RELATIVISTIC DESER-TYPE FORMULAE WITH LOWEST-ORDER 

CORRECTIONS 

In the lowest-order calculations only the first two terms in the r.h.s. of Eq. ( p3D con- 
tribute. Further, one can assume ipc{0) = / d'^k / {27r)^ ipc{M* ; k) = m~^/^0o in this approx- 
imation. Then, taking the real and imaginary part of Eq. (^3j), we arrive at 

Further, we can write 

Re (zTia) = 167r %+^-^^+^- {Ami; 0. ) 

Im(zri2) = -16vr (^)'^'(l - | W^.o.o(4m2; 0, go )p, (25) 

where T{s]p,q) denote the (dimensionless) S-wave rni scattering amplitudes and go is the 
relative momentum of the vr^vr^ pair at the threshold s = 4m^, with the magnitude given by 
the relation = m^o + Qq^- 
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We would like to emphasize that the second relation in Eq. (pSj) differs from an analogous 
relation given in Ref. though the starting equations (24) in the both papers coincide. 
In Ref. the magnitude of the three-momentum both for vr+vr" and vr^vr'^ pairs was set 
equal to 0. Consequently, neutral pions in the final state turned out to be off shell. Contrary 
to Ref. |22[ we deduce from the Cutcosky rule that neutral pions in Eq. (|25| ) are on-shell. 
Note that this discrepancy with Ref. p2| leads to different predictions for the corrections 
e.g., due to the mass difference m^± — 171^^0 in one loop (see below). 

Substituting Eq. (pSj) into (p^), we reproduce the lowest-order relativistic Deser-type 
formulae for the energy level shift and lifetime of the pionium 



AE 



r 

o TT+TT 

mi 



167r /2Am^\i/2 



mi \ m^ 



1 - 



2m7r 



IX- 



(26) 



A. Relativistic correction to the pionium w.f. 

The correction in the pionium lifetime coming from this effect is contained in the second 
term of Eq. ([23|). Namely, up to 0(a) terms 

^c(O) = / 7^ MM'; k) = 4^(1 - Coa), ^ = 0.381 • • • (27) 

(details can be found in Appendix A). Thus, up to order 0{a) in the local approximation 
the second term in Eq. ( [23| ) reads as m~^ T12 (p^ {1 — 2Coa). The term proportional to 
C^a in this expression induces the corresponding correction in the pionuim lifetime. Note 
that the value of this correction is determined by the expression of the unperturbed solution 
ipc (H) and hence depends on the particular choice of the instantaneous Coulomb part of 
the potential (for our case the Barbieri-Remiddi prescription (^)). We shall see, however, 
that in the full expression for the correction to the atom lifetime the term proportional to 
CqU disappears, indicating that the final results do not depend on a particular choice of the 
zeroth-order kernel. 



B. Correction due to the displacement of the bound-state pole by strong interactions 



This correction is induced by third and fourth terms in Eq. 



The calculation of this 



sort of integrals is carried out in a straightforward way was considered in Ref. Below 
we give the result of these calculations 



<Vc|(G'o')'G'oTi2|^c7> =Ti2^c(0) / 



{2n) 



^MM'';p){G^\M^;p)yGo{M*;p) 



o? mi 



+ ■ ■ 



(28) 
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< ^cliGori^c > = I -^MM'-^P)iGo\M*;P))"MM*;p) = 4 + " " " (29) 



<ijc\Ti2G,{G,y\iljc> =^c(0)Ti2 I -^^G,{M^-p){G,\M^-p))"ijc{M*-p) 



" (2vr)^- 

= -4 ^ + ■■■=< ^c|(Go ^)"GoTi2|^c > (30) 

and ellipses stand for higher-order terms in a. From Eq. (^4]) it is easy to see that the real 
and imaginary parts of integral ( pO| ) are down by small factors ^E^"^"* /m.,^ and r*^^Y^7r as 
compared to integral (pQ]). Therefore further we shall neglect (^) in our calculations. 

Note that Eq. contains the second derivative of the inverse free Green function Gq 
with respect to the bound-state mass. Hence this is a true relativistic correction arising 
from the BS treatment of the bound-state problem since in the nonrelativistic case the free 
inverse Green function is linear in the bound-state energy (cf. Ref. |P4|| ) 

C. Correction due to the exchange of Coulomb photons 



This correction stems from the fifth term of Eq. (^5]). The calculation of the correspond- 
ing integrals is considered in Appendix A. Below we give the results 

<,ljc\{G,')'6GT,,\^c> = TuMO) J -0^,-0y,MM''-,p){Go\M*-p)y6G{p,q) 



+ (31) 

a"^ mi 



< ^c\T^2SGT,,\i;c > = {Tuf (MO))' J 7^7^'^^(P'^) 



(27r)4 (27r) 

(Ina - 2.694) (T12)' + (32) 



where the term nonanalytic in the fine structure constant (containing Ina in Eq. (0)) comes 
from the infrared-singular one-photon exchange piece in the Coulomb Green function (Eqs. 

(0), (ii))- 

With the calculated integrals and the lowest-order relations (p4|), it is a simple algebraic 
task to derive, from Eqs. (^) and (^), the first-order correction to the pionium decay 
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r(^)=r«{(-^^) + (-2^0) + (1/2 + 2.694 -Ina)^^''' 



, El 

^ ' relativistic w.f. v 

strong Coulomb photon exchanges 



-(M*r«) 'Re<^c7|(l+Ti2Go)V^3(l + G'oTi2)|^c>} (33) 
where Ei stands for the energy of the unperturbed ground-state level Ei = —1/4 m.^a'^. 
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Comparing the Eq. (|33D with the corresponding expression given in Ref. p2[ (referred 
to as "second-order strong correction"), it is easy to see that our expression for the contri- 
bution of Coulomb photon exchanges contains an additional Ina term. The origin for this 



disagreement can be easily established. Namely, the authors of Ref. p2[ include an additional 
contribution from the so-called "constraint diagram" into this term. This diagram cancels 
explicitly the one-Coulomb photon exchange term, leaving only multiphoton exchanges. We 
have checked that, having merely discarded this term, after our calculations we come to the 



result numerically very close to that given in Ref. |^2[- The true result, however, cannot 
depend on the formalism used for the description of a bound state, either the BS equations 
or equations of the 3D constraint theory; So, one can ask whether the contribution regarded 
as the counterpart of the "constraint diagram" exists in the BS framework for bound states. 
The answer is yes, this diagram is contained in the "electromagnetic kernel" V3 in Eq. (p3D. 
The reason why we include this diagram into V3 rather than in the "second-order correction" 
is simple: this diagram is accompanied by the diagram of virtual photon exchange (see Fig. 
2d) neglected in Ref. The latter diagram also produces the Ina term which exactly 

cancels the corresponding term from the "constraint diagram" (below we shall discuss this 
in more detail). The remainder is regular in the fine structure constant. Consequently, we 
find more natural at the present stage to omit both these contributions on the equal footing, 
rather than to retain only one of them, namely the "constraint diagram". 

We would like to mention here that the sign and magnitude of the nonanalytic term 
appears to be exactly the same as in the nonrelativistic treatment of the pionium, indicating 
that, as one could expect from the beginning, the "electromagnetic kernel" produces the 
corrections which are analytic in a (at least in the lowest order). In the scattering theory, 
when the electromagnetic corrections are taken into account, the expression for the decay 
width corresponding to the second relation from Eq. (^) ImiTi2 is replaced by ImiTcc, 
where Tec denotes the scattering amplitude of charged particles at threshold in the presence 



of the Coulomb potential [0,^. Accordingly, it leads to the replacement of Imahad by 
Im Occ, where atad and acc denote the "hadronic" and exact scattering lengths of charged 
particles. However, assuming that the hadronic potential has a finite range denoted by R, 



the following relation between a^ad and a^c can be established |T3|J28|] 



1 1 2 , f2R\ If.. ,2R\ 
— = m — H series m powers of — ( 34 ) 

acc CLhad Tb KteJ TbK Tb J 

where is the Bohr radius of the pionium which is inverse proportional to a. So, up to 
the logarithmic terms, 

/ 4Re tthad , \ / AE(^) , \ , , 

Im acc = Im ahad (^1 Ina j = Im ahad [I ^ — ma j (35) 

where in the second relation we have used the Deser formula for the atom energy level dis- 
placement. Consequently, the decay width is modified according to F ^ T (1 — AE^^Hna / Ei) 
(cf. Eq. (ID). 

In addition, we would like to note that the logarithmic term which, just as in our ap- 
proach, is of the second order in strong interactions (~ (T'12)^) was found by Roig and Swift 
p9| . The authors of Ref. ||2^ have studied the electromagnetic radiative corrections to the 
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TTTT scattering and discovered the term proportional to Inp in the amphtude. When sub- 
stituted into bound-state integrals, after rescaling, as usual, the integration momenta by 
P ~^ IVi this term produces the Ina correction, as in Eq. (p3|). 

So far our treatment of the pionium decay width has been incomplete. Now we turn to 
the calculation of corrections induced by the last term in Eq. ( p3D containing the "electro- 
magnetic" kernel V3. 



D. Mass shift and radiative corrections 

Into the kernel V3 we include the following diagrams: the diagram with the "residual" 
photon exchange (i.e. the virtual photon exchange minus Coulomb potential). Fig. 2, the 
self-energy corrections to outgoing pion legs. Fig. 3, the vacuum polarization diagram. Fig. 
4 and the vertex corrections. Fig. 5. The contributions containing low-energy constants and 
tadpole terms are included into T12. Thus, we take here the advantage of the arbitrariness 
in splitting the potential and include all terms into T12 which are a priori known to have a 
smooth momentum dependence on the bound-state scale 7. Only the potentially "danger- 
ous" terms which contain the photon propagator with vanishing mass are to be treated with 
the bound-state equation. 

In this section we are concerned with the first two terms in V3. They read as 

1. Residual photon exchange (Fig. 2) 

V^-Vc = te\P + P + q)^{P -p- q)uD^''{p - q) - Vc{p, q) (36) 

where D^" denotes the photon propagator. The calculations are most easily carried in the 
Coulomb gauge with 

D'%k) = -j-^, D'\k) = D^\k) = 0, DHk) = -[S^^-j^)^^ (37) 

The corresponding matrix element equals 

Re < ^c\{V^ - Vc) + TuGoiV^ - Vc) + {V^ - Vc)GoTu + TuGo{V^ - Vc)GoTu\^c > 

It is easy to observe that Re < ipcliV'y — Vc)\'^c >= at the bound state energy. Conse- 
quently, the first term in the matrix element vanishes. The following two terms are equal 
for symmetry considerations. Further, we completely neglect the fourth term in this matrix 
element (Fig. 2d). This term, as we have mentioned before, contains the contribution of 
"constraint diagram" given by Vc- However, it is obvious that the D^^ component in the 
expression of V-y which has exactly the same infrared singular behavior as the Coulomb po- 
tential leads to the same Ina nonanalytic term in the lifetime. In the expression of V3 these 
nonanalytic terms cancel and the remainder is analytic in a (at least in the lowest order). 
For this reason we find more safe to neglect the combination — Vc rather than Vy alone, 
as in Ref. 



The remaining term in the expression of the matrix element is exactly of the form of 
the "retardation correction" discussed in Ref. [^. To demonstrate this, we note that 
ISi/j^ >= GoiVy — Vc)\4'c > gives the first-order perturbative correction to the bound-state 
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w.f. due to the retardation effect (i.e. the difference between Vy and Vc). The correction 
induced in the pionium hfetime due to this effect is given by 



r 1- 



1 



^Re[2^c(0)Ti2 5^^(0)] 



r 1 + 



(38) 



where the second relation was derived with the use of Eq. (|2^. From Eq. (^) it is apparent 
that the kernel — Vc- is responsible for the retardation correction in the pionium lifetime. 
Finally, the matrix element we are looking for can be written in the following form 



= Re 



2 X ^c(O) X X 



(39) 



where 



A. 



-26 



Go{M*- p){P + p + q)^{P -p-q)u D^^^ip - q) 



A, 



c 



d'^p 
(2^ 



G'o(M^p) Vc{p,q) 



(40) 



2. The kernel corresponding to the insertion of self-energy graphs into the outgoing pion 
legs (Fig. 3) 



where 



n(/) 



d^k 



+ 



n(f -g) 



n(f + 9) 

T + qf -ml ' (f - qY - ml 



(41) 



(27r)4 {l-kf-ml 



{21 - k)^{2l - k), D^'ik) 



ti{l)=Ii{l)-Yi{f 



Z{1) 



[2 _ ^2 ' 

The corresponding matrix element can be rewritten as 

-Me =< ^cl Tu {Z+ + Z_) (1 + G^T^^) > 



Z^{l) = Z{-±l) 



(42) 



(43) 



We neglect here the term which is of second order in the strong interaction amplitude. Then, 
the sum of matrix elements reads as 



+ Mt. = Re 



2 ^c(O) Ti2 



(^(A, -Ac- jZ^q) - jZ_{q))MM*; q) 



(44) 



We would like to emphasize here that the coefficient | which emerges naturally in front 
of the self-energy term in the r.h.s. of Eq. (^) ensures that the correction to the decay 
width calculated from (H3) is gauge invariant. Moreover, it is well known that in this 
particular combination of vertex and self-energy diagrams the "photon mass" disappears in 
the calculated width without inclusion of the "soft photon emission" terms, which seem to 
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be rather awkward in the context of the bound-state problem. In its turn, the origin of 
the coefficient | can be traced back to Eq. (^), where, as was mentioned before, the kernel 
V contains only a half of all possible self-energy graphs attached to the external pion legs. 
Note also than in the S'-matrix elements the origin of emerging the coefficient | is quite 
different. In the latter case one takes into account all self-energy graphs in the external legs. 
However, restricting these matrix elements on mass shell one encounters an expression of 
the type 0/0 which can be tackled, introducing an explicit "smearing function" in the initial 
Lagrangian. Passing then to the limit when the "smearing function" tends to the unity, one 
discovers the coefficient | which multiplies the contribution from the self-energy diagram 
(for a detailed discussion see e.g., Consequently, though the coefficient | emerges in 

the bound state and the scattering problems from different sources, one ends up with the 
same expression in these two cases, being gauge invariant and infrared finite at threshold. 



The integral in the r.h.s. of Eq. (^4 ) can be easily evaluated in the Coulomb gauge, 
bearing in mind that the final result is gauge invariant. 

In the calculation of the contribution from one can use the fact that the D^^ component 
of the photon propagator in the Coulomb gauge contributes only in order 0{aHna) in the 
decay width and thus can be neglected. This considerably simplifies the calculations. The 
result reads as (details can be found in Appendix B) 

^A,^c(M^g) (45) 
(2^(2^0^^°^^ ^'^""^ ''^ ''^"^'^^ 



where the dimensional regularization was used to handle the ultraviolet divergences, and 



with n being the dimension of space and /x the mass scale used in dimensional regularization. 
The contribution containing Ac can be trivially carried out 

"^'^ KcMM*;q) (46) 



(27r) 

A ^ N4«m^e2(w(p)u;(g))i/2 ^ ^,^ATiam^(po 

(27r)^ (27r)* \p — q) q^^Y 



d^p d^q 111 47ram^0o 0o 

(27r)2 (27r)2 w(p ) V2 p 2 + y 2 + ^2)2 = ^ 



The pion self-energy graph (^2]) calculated within the dimensional regularization scheme 
in the Coulomb gauge is given by 
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where A = — ml + iO and 



2 + ^lnr^(0 L = 0(/-) at small ir 



l\ \w{l) + \l 



1 /.I 



JO 



/(/^A)= / dx I du ^ (48) 

x(l — M^)/ ^ + xm^ — (1 — x)A 



Note that the terms in square brackets in the r.h.s. of Eq. ( ^7D are of higher order in 
\l I and/or A as compared to the second term in the round brackets. Consequently, in the 
calculations in the lowest order in a the term in square brackets can be dropped, since it 
contributes only in order a^lna. Then we immediately obtain 

and 

= Re(2 X {MO)r X T,, x {-e')^-^^N, - (50) 
Putting things together, we finally obtain 



3a ,^ ^ ct 



M^ + Mj: = 2ReTi2 ^ [—N, + Coa- — ) (51) 

V47r ZTT/ 

which induces the corresponding first-order correction in the tt+tt^ atom decay width 

r^r(|iv. + 2C„«-^) (52) 

The above expression is of course ultraviolet-divergent. It is well known that, along with 
the diagrams contributing to this expression, one should consider the four-pion Lagrangians 
containing (divergent) low-energy constants in order to cure this ultraviolet divergence. This 
will be done below. We would like to mention here that the term 2Coa from this expression 
cancels with a similar term coming from the atom w.f. in Eq. (P^D, and the final result for 
the decay width does not depend on the initial approximation chosen for the Coulomb w.f. 
of an atom (as it should be). Further, the term —a/n exactly coincides with the result given 
in Ref. obtained from the same set of diagrams in an arbitrary covariant gauge. This 



provides an independent check of the gauge invariance of our result also for noncovariant 
gauges (and, in particular, for the Coulomb gauge). 

Below we would like to discuss briefly the connection of our result with the "retardation 



correction" given in Ref. pO[. As we have mentioned above, the matrix element Ai-y gives 
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exactly what is called "retardation correction". Note, however, that our result differs some 
what from that of Ref. 
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Namely, in this paper the virtual photon exchange diagram 
corresponds to the Wick-Cutcosky model, whereas we have used the pion-photon vertex 
which emerges in scalar electrodynamics. The ultraviolet divergence which occurs in our 
result is a consequence of the choice of the pion-photon vertex. Thus, strictly speaking, the 



present results, and the results of Ref. refer to different physical models, and cannot 
be directly compared. Further, as we have seen, the contribution from Ai^ alone is gauge- 
dependent and should be combined with the self-energy diagrams to yield a gauge-invariant 
result. Moreover, in gauges other than the Coulomb gauge, individual contributions from 
Aiy and A^s contain a nonanalytic Ina dependence which cancels in the sum. Reference 



20| which mimicks the Feynman gauge calculations, does not contain a nonanalytic term. 

Let us now evaluate the contributions from local four-pion Lagrangians. As we have 
mentioned before, it is appropriate to include these terms, which are smooth functions of 
external momenta, into the definition of T12. According to this convention, the transition 
amplitude %.+^~^^o^o from Eq. (^) can be written as %.+^~^^o^o = + 7^, where Ti 
denotes the isotopically symmetric "strong" vrvr scattering amplitude with the mass of the 
isotriplet taken equal to the charged pion mass, and T2 includes the effect of isospin breaking 
as well as the terms with low-energy constants from the four-pion Lagrangians. From Eqs. 
and (|52D we come to the expression 



r(2) 



Ztt 71 Ti 



(53) 



which displays only the electromagnetic and mass shift corrections 
As we noted before, in Ref. 



2^ the amplitude T2 was evaluated at the off-mass-shell 

the amplitude emerging here is 



point for vr mesons. However, as we see from Eq. 
restricted on mass shell for all external particles, and we use this prescription hereafter. 
Moreover, the explicit expression for this amplitude calculated within ChPT has become 
recently available [|31[, and in the following we can use the expression given in Ref. as 
granted. All what we have to do is to extract from the amplitude of Ref. the terms 



which we have already taken into account through the bound-state equation (vertex and 
self-energy corrections, i.e. only the ones which are taken into account in the model of Roig 



and Swift |2|) 



The calculations in Ref. ||3T]] were carried out in the Feynman gauge. However, as we 



mentioned before, the combination A^ — \Zj^ — \Z_ we are concerned with is gauge- 
invariant, and we can safely use the results of Ref. Thus, we can identify (see Eqs. 
(4.9)-(4.11) and (4.19) from Ref. [§T|) 



- 2e^{s - 2ml)G^ 



where 31 



G+_-y(s) = —i 



is) - e^A 



A. 



3e^ 
167r2 



87r2 I A2 



47r2 



1 



(27r)4 (g2 _ A2)(g2 _ 2g • p+)(g2 + 2g ■ p 
1 



V± 



P 
2 



(27r)"(g2_m2)((g_/)2_^2 



2V 



(54) 
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and we have introduced the photon "mass" A to regularize the infrared-divergent integrals. 

2 2 

The amplitude 7^+^-^^o^o can be easily found, subtracting — — defined 



by Eq. ( [5^ ) from the total amplitude given in Ref. 



- 327rr^ 



+ 7r — ^TT^TT*^ 



s — mzo 



BR{s,t,u) - CR{s,t,u) 



with 31 



(55) 



BR{s,t,u) 



s — mz.0 



+ 



+ 



1 



12F4 
1 

12F4 



A2\ 2 



X{t,ml,mlo) 



Al + t{s-u) 



A, 



2 \ 2 



3(M-2m^ + ^) + 



X{u,ml,mlo) / . 2 



J+o{u) (56) 



Cii{s,t,u) 



327r2F4 



4A, 



{Ami - 7rniomi + 5mlo) + e'F\-6 -QN, + IC 



+ 



+ 



327r2F4 

247r2F4 
1 

1 

487r2F4 



-f ^ - ^(2^2 - m^o) + m^o/s + e^F^/Cf 



A. 



967r4F4 



1 1 

t u 



A. 



/i(s-2m2o)(s-2m^) + 



i (11,2 _ ^2 _ ^2) ^ ^ / 2 _ ^^2^ ^2 ^ ^2^^ ^^2 ^ ^2) 



487r2F4 



r2[(t-s,)2 + (M-so 



Here 



±0 



jr±o 



4Z\- 
3 + -)*. 



32Z 
~9~ 



8Zk2 + 3k3 + AZki - 2(1 + %Z)kQ - (1 - 8Z)A;8 



2 I 2 



A. 



2 2 

TV TV ' ' 



In 



A. 



2e2F2 



(57) 



(58) 



and li, ki denote the finite, renormalization scale independent low-energy constants |^ 



Note that in the expression for Cr the term containing the photon mass has been explicitly 
cancelled (cf. Ref. ^Tj]). The term proportional to A^^^ in this expression exactly cancels the 
ultraviolet divergence which appears in Eq. (^2]). To demonstrate this, we note that in the 
presence of this term only the lowest-order scattering amplitude is modified as (see Eq. (|57[) 



s - mlo ( ^ 3a; 



p2 



p2 



1 

An 



(59) 
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The modification in the vr+vr" atom decay width is twice as large, and this cancels the term 
in Eq. (^). Consequently, we can merely discard the ultraviolet-divergent quantities 
simultaneously in Eqs. ( ^2] ) and (|57|). 

It is convenient to expand the rest of the vrvr scattering amplitude, which is ultraviolet 
and infrared stable, in powers of A^^ near threshold. Below we present the result of this 
expansion retaining only the terms proportional to A^^ and e^F^ which are sufficient up to 
an accuracy required. Denoting the corresponding amplitude by T+^^-^^o^o, we find 



- 327rr(? _ 

TT'TT 



7r~7r — ^TT'TV^ 



'iml 



ml 



+ A, 



F2 487r2Fi 



F2 327r2F4 



87 16- - 
11 + -h + y/2 - /3 



12L 



32772^2 



(60) 



where the charged pion decay constant F^^ is related to the parameter F entering into the 
Lagrangian through [pl|] 



^'0 
167r2F2 



(61) 



The low-energy constants U and fcj in this equation are fixed on the renormalization scale 
/i2 = m^o according to 



Vi 



327r2 



k + \n 



A:[(^) 



327r2 



ki + In 



(62) 



where rji and cTj are constants. However, to make the comparison with the calculations 
carried out in the isotopically symmetric case, it is necessary to bring the normalization 
scale to m2. This induces the change in the second term of Eq. ( |60D 



487r2Fi 



;i + + 3r3- 12/^4) 



487r2Fi 



- - - 1Q?7?2 

1 + 4/1+3/3-12/4)+ " 



327r2Fi 



(63) 



After this rescaling the first term gives the isotopically symmetric "strong" amplitude 
Ti and the remaining part corresponds to T2. Using then Eqs. (|53|) and (l^); it is easy to 
"read off" the first-order mass shift and radiative corrections in the h^tx' atom decay width 



2A. 

3m2 



1 + 7^(1 + 4^1 + 3^3-12/4) + ;^^ 



m: 



967r2F2 



h + 12/4 



(64) 



= 2Coa - - + T^(-18 + 3/Cf - /Cf ) (65) 

It is interesting to note that the Deser-type formula with account of the first-order mass 
shift and radiative corrections can be rewritten in a simple and transparent way. Namely, it 
is well known that the scattering amplitude of charged particles develops a pole at threshold. 
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which corresponds to the long-range Coulomb interactions in the initial state. Thus, at the 
threshold we can write fST 



16 q 

where q is the c.m.s. relative three-momentum of charged pions. 

Then, the following simple expression, valid in the lowest-order approximation, is ob- 
tained for the TT+vr" atom decay width 



GATTml V niTT J \ 2m. 



{ReA 



thr J 



X (^1 + (^-- -^j + (-2Coa) + (+2Co«) + (1/2 + 2.694 - Ina) + 5 j (67) 

where the electromagnetic and mass shift corrections, being excluded from 6, are completely 
taken into account in ReA^~'°°. Thus, the correction factor 6 displayed in Eq. (p7|), includes 
the contributions from other sources, e.g., vacuum polarization, finite size corrections, etc. 



Further, equation (|67|) demonstrates explicitly the cancellation of 2Co« terms which depend 



on a particular choice of the initial approximation for the relativistic Coulomb w.f. 

It is worth noting that the quantity Rey4^~'°° is not proportional to the conventionally 
defined vrvr scattering lengths which acquire an additional finite contribution due to the 
emission of real soft photons (see Eq. (5.17) from Ref. ||3l]]). Thus Eq. ( |67D demonstrates. 



that in the presence of long-range Coulomb force the pole-removed real part of the scattering 
amplitude rather than the scattering length enters into the expression of the first-order 
corrected Deser-type formula for the decay width. 

Below we shall briefiy discuss the comparison of our results with those obtained in Refs. 



2^ , |21[| . As we have noted before, the main difference between our work and paper p2[, where 
the corrections to the pionium decay width are also evaluated in ChPT, consists in the fact 
that we argue the necessity of a different kinematic prescription in the calculation of the vrvr 
scattering amplitude entering into the Deser-type formula. It is a completely on-mass-shell 
amplitude which naturally emerges in our calculations with neutral pions having small, but 
nonzero relative momentum go in the final state. By contrast, in Ref. both charged and 
neutral pions have zero relative momenta, and, consequently, neutral pions in the final state 
are slightly off-shell. Different kinematic prescriptions lead to different predictions for the 
pionium decay rate in Ref. ||22| and the in present work. Further, in Ref. I^T] the radiative 



corrections to the tttt scattering amplitude were evaluated in the Roig and Swift model, with 
an explicit cutoff. The cutoff parameter was chosen to be equal to the p-meson mass. The 
author also presents the calculations carried out in the p-meson dominance model where 
the integrals are ultraviolet convergent and the explicit cutoff is not needed. Thus, a direct 



comparison of the results of Ref. |^ with the calculations carried out on the basis of ChPT 
is not possible. Note, however, that in Ref. the contribution from real photon radiation 
is also included into the scattering amplitude. 

The mass shift and electromagnetic corrections have been evaluated in the framework 
of the nonrelativistic scattering theory approach |[T^JT5[|. Strong vrvr interactions in this 



approach are described by energy-independent local potentials. It turns out that the effect 
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of the mass splitting in the pionium hfetime is opposite in sign as compared to the case 
of field-theoretical calculations. It is obvious that this sign depends on the choice of the 
"reference mass" corresponding to the case of the isotopically symmetric world. In particular, 
if the reference mass is chosen equal to the charged pion mass [0, then this effect turns 
out to be negative (~ —3%). If one chooses the neutral pion mass to be the reference 

+3%). 



mass 



Tj], then this effect changes its sign 



our calculations based on the chiral Lagrangian. 



We observe the opposite situation in 
Thus one can conclude that the energy- 



and mass- independent local potentials used in Refs [p!4| , p!5| might not provide an adequate 
description of the vr+vr" system in the nonrelativistic limit in dealing with the sophisticated 
issue of the isospin breaking effect in tttc interactions. 



E. Correction due to vacuum polarization 

The vacuum polarization due to the virtual electron-positron pair contributes in order 
to the pionium decay width. However, this effect is amplified since a small electron mass 
rrie is present in the denominator. 

In the instantaneous approximation the photon propagator is modified by the vacuum 



polarization effect as follows [18 



1 1 a , ^ /V 

where 

J Ami S + S\ S / \ S / 

The perturbation potential which is responsible for the vacuum polarization effect is given 
by 

16i 

V^ac = -^mla^ hac{-{p~ q?) (70) 
and the corresponding matrix element from Eq. (|33|) is equal to 

M^ac = Re < MVvac + TisGoKac + Kac^oTis + TisGoKacGoTis | V^C > (71) 

The first term in this expression vanishes at the bound-state energy which is below the 
elastic threshold. As in the calculation of electromagnetic radiative corrections, we neglect 
the fourth term (Fig. 4d). Thus, the matrix element can be written as follows 



Mvac = Re 



"^9 r d^n d^n 

mla' zTi2 MO) / 7^7^Go(M^; p)4,,(-(p - g )2)^c(M*; g) 



3 " - -^ov y y (27r)4(27r 



(72) 



With the use of Eq. (pOD and after integrating over relative energy variables the integral in 
Eq. (|72D takes the form 



coo 



I d ( ) [ -inam^cpQ 

Ami ^^^""^ J i27iyi2n)^ wip)iwiq)y/^{p^ + j^){q^ + -f^yis + {p-qy) ^ ^ 
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In the calculation of the integral over d^q we can replace the smooth factor (u'(g'))^/^ by 
its value at g =0. Then 



(27r)3 («;(g))V2 (^-2 + ^2)2 (, + g^^^i/2 ^-2 + + 



+ 



(74) 



The remaining integral over d^p can be computed analytically. Rescaling the integration 
variable s in Eq. (|73D , we finally arrive at 



M 



vac = ^77"^ — ReTi2 60 

lb nip 



(75) 



where 



roo dsis - ly/^ 1 

1 + Yg )f^^c{s) X 



dsis-iy/^ 



1 



1 -1 



0.6865- 



(76) 



2m, 



7r(27 + yi) 
+ 6'(a/s - + 7) 

7 



6'(m^ - 7 - ^) 



7 + 



(m2 _(^ + yi)2)i/s 



arctg 



(^^-(7 + v^)^)^/^ 
7 + \/i 



7 + 



In 



7 + + ((7 + V^)' - ^71^) 



2\l/2 



{ml — 7^)1/2 



arctg 



((7 + v^y -miy/^ V 7 + - ((7 + v^y - miy/^ 

(^2 _ ^2)1/2. 



7 



(77) 



and = m,r/(2me), 7 = 7/(2me). 

Using the relativistic Deser formula in the lowest-order approximation (pi]) , it is easy to 
observe that Eq. (^) leads to the following modification in the vr+vr" atom decay width 
due to the vacuum polarization effect 



(78) 



Note that in Ref. 18 



the vacuum polarization correction to the pionium lifetime was 
calculated only with account of discrete spectrum transitions. Thus, our result is a general- 
ization of that from Ref. 



m 



The NRQED based calculation of the vacuum polarization effect in the pionium lifetime 
has become available recently p3[. We find that the analytic expression of the so-called 
0-Coulomb term in Ref. coincides with our result up to the relativistic kinematic factor 
in the w.f. whose presence is due to the choice of Barbieri-Remiddi kernel. Numerically the 
effect of this factor, which contributes in higher orders in a, is very small. In Ref. pSf the 
results for 1-Coulomb and multi-Coulomb contributions are also given. These contributions, 
which formally are of higher order in a, would emerge in our calculations as second-order 
perturbative corrections to the pionium lifetime. 
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F. Finite size correction 



In the presence of a pion loop the bare VPP vertex is modified. According to Ref. 

Mt) = 1 + - 4m^) J+_(t) + ^^(^6 - ^) = 1 + ^Mt) (79) 

In the instantaneous approximation the perturbation potential is given by (see Fig. 5) 

V, = S^e^nl '"'fy-f^ (80) 
[p-qy 

and the corresponding matrix element reads as 

Mf =< ^Pc\Vf + TuGoVf + VfGoTu + TuGoVfGoTu\^c > (81) 

Again, the first term in this matrix element vanishes at the bound-state energy, and we 
neglect the fourth term. The remainder is then given by 



Mf = Re 



2V^c(0)Ti, /(^(^Go(M;p) ^c(M;g) 



(82) 



The integral over d^p diverges in the ultraviolet region. This stems from the fact that the 
diagram depicted in Fig. 6a, is ultraviolet-divergent as a whole, though the subdivergence in 
the VPP vertex has already been removed by an appropriate counterterm depicted in Fig. 
6b (this counterterm is implicit in Eqs. ([79|) and (^2])). Thus, a higher-order counterterm 
shown in Fig. 6c is needed to cancel the overall divergence in Fig. 6a (and, consequently, in 
Eq. (p2D ) . It is obvious that this divergence is removed analogously to that from subsection 
D (Mass shift and radiative corrections), and we shall not further dwell upon this question. 

To simplify the calculations in the relevant low-t region, instead of Eq. (|79D we use the 
well-known monopole parameterization 

bFy{t) = ^'v^ = \< >v (83) 

with the same < >y as in Eq. ([791). The integral in Eq. (p^ is then convergent and can 
be easily evaluated, resulting in 

M (27r)3(27r)3 W(p)(w(g-'))l/2(p 2 + y)(^2 + ^2)2(^2^ + (p'_^)2) ^ ) 

With the use of Eq. (|74D and the inequalities nv » "^tt >> 7 the integration of Eq. (|84D 
gives 

_^2 2^^ ( d^P \ 00 e^m^ 4^ 

"""^my^i {27,fw{p){p^ + l^){P^ + l^l) my^27rVy K " ' ^^^^ 

and 

>,, = -^ReT„i^^!!|lu^4 + ... (86) 

The modification of the tt+tt^ atom decay width due to the finite size effect is given by 

r^rfi-l^^!^inl4) (87) 
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IV. NUMERICAL RESULTS AND DISCUSSION 



In this section we present the numerical results on the lowest-order corrections to the 
pionium decay width. To this end we combine various corrections obtained in the previous 
section. As we have seen, the correction due to the relativistic modification of Coulomb w.f. 
cancels with the corresponding piece in electromagnetic radiative corrections. In the final 
result we cancel these corrections explicitly. 

Below we give a list of analytic results on the lowest-order corrections to the tt+tt^ atom 
decay width 



167r /2Am^\i/2/. Am^\i/2 2^2 j2 m . 

[al - ao) 00 [1 + Os + Oc + Om + Oem + Ovac + Op) 



9 V / V 2m 



(88) 



where and denote the tttt scattering lengths in the isospin-symmetric case, with the 
charged pion mass taken to be the common mass of the pion isotriplet. 
5s is the correction due to the displacement of the bound-state pole by strong interactions 
(see Eq. (0)) 

9 Ai?(^) 

= —ET- = -5-47 X 10-3 m^(2a° + al) (89) 

O -C/1 



5c corresponds to the correction due to the Coulomb photon exchanges (Eq. (pSD) 

5c = (1/2 + 2.694 - Ina)— ^ = 3.95 x 10"^ m^(2a° + al) (90) 

El 



5m stands for the correction due to the m^± — m^o mass difference (Eq. (|5iD) and 5ern 
corresponds to the electromagnetic corrections without 2Coa term (Eq. (^). The quantity 
5yac denotes the correction due to the vacuum polarization effect (Eq. ([78|) ) 



5.ac = 4«'—x 0-6865 (91) 
lb rrie 

and 5f corresponds to the finite size correction (Eq. (|87D) 

5p = ^ml<r^>-y ln(lm^<r^>^) (92) 

To make numerical estimation of the above-listed corrections, one has to substitute the 
values of low-energy constants into these expressions. For the constants U we take the 
numerical values from Ref. h = -2.3 ± 2.7, k = 6.0 ± 1.3, I3 = 2.9 ± 2.4, h = 4.3 ± 0.9. 



The constants ki are more difficult to estimate. In our paper we use the values from Ref. |31 
based on a rough estimate at the scale coinciding with the p- meson mass |/c[(mp)| < 

This estimate yields ^ /Cf ° = 1.8 ± 0.9, ^ /Cf ° = 0.5 ± 2.2 0. Large error bars in the 

low-energy constants /Cf°, /C^°, in turn, do not allow one to calculate the electromagnetic 
radiative correction to the atom decay width with a high accuracy. Other input parameters 
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in our calculations are the S-wave tttt scattering lengths: «[] — 0.217m~^, Oq = — 0.041m~^ 
calculated in ChPT and the e.m. charge radius of pion < >v= 0.439 Fm 

Substituting the above values of the input parameters into the expressions for various 
corrections to the decay width we obtain the results collected in Table I. As we observe from 
this table, the largest correction in the decay width is caused by the mass splitting effect in 
accordance with the result of Refs. |ll4| , p^p2| . Our result for this effect has the same sign as 
the result of Ref . |^ , but is different in magnitude. The reason for this difference is traced 



back to different kinematic prescriptions used for the tttt scattering amplitude in the Deser 
formula (see the discussion in the text). The sign of the mass splitting effect obtained in the 
nonrelativistic scattering theory approach turns out to be opposite as compared to 



our result, and is of the same order of magnitude. 

Our last remark concerns the effect of the rud — ruu mass splitting on the pionium decay 
width. It is well known that in the one loop order this leads only to the shift in the neutral 
pion mass [^. Since in our calculations of the on-shell tttt scattering amplitude we use the 
physical values of the pion masses, with the mass difference caused both by — 7^ 
and electromagnetic corrections, the resulting mass splitting correction includes both these 
effects. [] 
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APPENDIX A: 

In this appendix, we present the calculation of various integrals which appear in the 
first-order correction terms. We start from the evaluation of the relativistic Coulomb w.f. 
at the origin (Eq. (|27|) ). By carrying out the po integration with the use of Cauchy theorem, 
this quantity can be written as 

f d^p p^ Airam^ \ 

^ ~ / TTTT^ , ,^...1^, 1/2 . ,^...1.., ' TZTT , .,2\2 \^^) 



In the lowest-order approximation in a we replace the factor ^ + 7^)^ in the denominator 
by (p and obtain 



^c(0) = ^(l-Co«) + --- (A2) 



We are indebted to Prof. J. Gasser for the clarifying discussions on this problem. 
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where 



9 roo 



dp 



Vr Jo (1 +p2)l/4(l + (1 +p2)l/4)(l + (1 +p2)l/2) 



0.381 



(A3) 



Next we turn to the calculation of the integral which is present in Eq. (|3lD . This integral 
contains the Green function SG{p, q) which corresponds to the exchanged Coulomb photon 
ladders and, according to the Eqs. ([T0|) , (pTf ), is given by 



6G{p,q) = i{w{p)w{q))^^'^ 

1 



mq)-sip)siq[ 



8 d 



M* DM* 



Go{M*;p)GoiM*;q) 



^{p,q) = IQ-Km.j.a 



R{p, q) = 25- 



{p- g)2 

1/2 



+ 



;S{p)S{q)R{p,q] 



(A4) 



[Sip) + sm + ■ 



where ellipses stand for the higher-order terms in a. 

Substituting this expression in the integral from Eq. (|3l|) and carrying out po, go inte- 
grations, we obtain 



^1 = / ,fjM M*-.P){Go\M*;p)y5Gip,q) 



(27r)4(27r)4 
d^p d^q AiramT^c, 



(A5) 



(27r)3(27r)3(p2 + 72)3[ 

where we have used 
dpo 



Hp,q)-S{p)S{q) 



(u;(g))V2g'2 + ^2 



(Go(M^ p)y{G^\M*; p)yGo{M*; p) 



3(M 



*\2 



In the calculation of 3D integrals containing the function g*), we use 
d^p I I 11 11 



+ 



(A6) 



and 



(27r)3 — g )2 2 + ^2^3 47ra;3m3 (g 2 + ^2^ 87ram7r (g 2 + 72)2 
d^p lR{p,q) f^dpf^ff ^/// 

-{Ji{p;q) - Ji {p;q)) 



where 

J[{p\q: 



d^p 



1 



1 



(27r)3 2 _|_ ^2-)3 g'^p _|_ ^^2q, 2(^p' 2 _|_ ^2-)(^g'2 _|_ ^2-)(^]^ _ p-)S 



(A7) 



(AS) 



(A9) 



The integration over d^p can be carried out with the use of Feynman parameterization. We 
obtain 
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3 1 r^dx{l-xf 3 1 f^dx{l-xy 

0\ \P\l) — o„„,fi™?i 1^1 , L ,3/2 ,5/2 



d± = \-x+^{\±pf (AlO) 



and 



J7'"( • / '^^^ ''^ ^ (All) 

Substituting Eqs. (|A11| ) and ( |A1CI| ) into (|A8| ) and carrying out the integration over dx and 
dp, we finally obtain 

- ' ' + „ ' ' (A12) 



(27r)^ (^2 + 7^)'^ 47ra;^m^ (^^ + 7^) Svram^ (g^ + 7^)^ 

With the use of Eqs. (|A^), ( [A7|) and ([A12| ) the integration in Eq. (|A5|) is trivially carried 
out, resulting in 

Substituting this result back in Eq. (|3lD , we readily obtain the final result given in this 
equation. 

Next we turn to the calculation of the integral which is present in Eq. (|32D 

where and correspond to the "nonderivative" and "derivative" terms in Eq. ( [A4D . 
Carrying out the integration over the relative energies with the use of Cauchy theorem, I2 
can be written as 

2~~16 7 (27r)3(27r)3 {w{p)w{q)y/^ (p ^ + 72)(g 2 + ^2) ^P^l) I ^) 

I2 receives the contribution from (a) one-photon exchange, (b) multiphoton exchanges con- 
centrated in Inip^q) and (c) the rest, proportional to the function R (see Eq. ( |A4D . Below 
we shall evaluate these contributions separately, 
(a) One Coulomb photon exchange 

, _ r d^p -iiram^ 1 1 _ iam^ f 3^ ^'^(r) 

~y (27r)3(«;(p)t/7(g))i/2 (p'2 + ^2)(^2 + y) (p'_^)2 - 4 J ^ ^ 

where 

'^^'■^ " / (2^ (z/;(p))i/2(p2 + ^2) (Al'^) 

Using exponential parameterization, the integration over d^p in Eq. (|A17| ) can be carried 
out, resulting in 



27 



4m 



— m(1 — a;)7^ — uxml 



^^^^ = Q 3/2r^l^ / duu~^^^ / (ia;a;~^/^ exp 
87r'^/^r(i) Jo Jo 

Substituting Eq. ( [A18| ) into Eq. ( [A 161 ) and integrating, we obtain 

^1/4(1 _^)-l/4^-3/4^^-3/4 



(A18) 



laniT, 



2a 



"327r3/2r2(i- 



(ir / (ixi / dx2 



T 



Vo " (72 + (m2 -72)(rxi + (l-r)a;2))i/2 



(A19) 



Note that one can not directly assume here 7 = in the denominator, since the integral 
over the Feynman parameters diverges in this limit. In order to overcome this difficulty, we 
split the integration area into two domains according to 

rl rl rl rl r-2 rl/p 

/ dxi dx2f{xi,X2)= pdp dtf{pt,p{l-t))+ pdp dtf{pt,p{l-t)) 
Jo Jo Jo Jo Ji Ji-i/p 

(A20) 

Performing explicitly the integration over dp in the first domain, we obtain 

^-1/4(1 _^)-l/4^-3/4^^_^)-3/4 



/2a(l) 



167r3/2r2(i)(m2 -72)1/2 7o 



1 /■! 

dt / dr- 



(rt + (l-r)(l-t))i/2 



xln [r\{ml - ^Y'irt + (1 - r)(l - t))'/' + (7^ + (m^ - l')irt + (1 - r)(l - t)))^^)] 



zalna ia 



IGtt IGtt 



21n2 + 



Cl 



+ 



27rV2r2(i/4) 

where ellipses stand for the higher order terms in a and 

.1 rl ^-1/4^1 _^)-l/4^-3/4^1_^)-3/4 



(A21) 



Cl 



dt / 

Jo 



(rt + (l-r)(l-t))V 



2 ln(rt + (l-r)(l-t)) = -40.374- •• 



The second integral converges when 7 — in the denominator, resulting in 

iac2 



/2a(2) 



4^3/2r2(i; 



(A22) 



(A23) 



where 



C2 = \ t dTT-^'Hl - t)-^I^ dp Inp (p - \)-^l^{r + (1 - r)(p - 1))-^/=^ = 0.288 ■ ■ ■ 
4 Jo Ji 

(A24) 

(b) Multiphoton exchanges 

In this contribution we can safely replace the smooth factor in the denominator 
{w{p)w{q)y^'^ — > m^r. Then 



d^p d^q Ir{p, q) 

(27r)3(27r)3(p2 + y)(g-'2 + ^2) 



—ma 



p 



(A25) 
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where according to Eq. ([lOD 



1 



X 



(27r)3(27r)3(p2 + ^2)(^2+^2) 

1 



(A26) 



.{p — <iYp + m~2a~2(p 2 + 7^)((f ^ + 7^)(1 — pY m~'^a~'^{p'^ + 7^)(g ^ + 7^) 
Introducing Feynman parameters and carrying out the momentum integration, we obtain 



167r2 



1 



dx 



2 Jo (d^/Vrff) 



(A27) 



and 



l2b 



IGtt 



(A28) 



(c) Factorizing integrals 

The integral containing the function R is evaluated in the straightforward manner, re- 
sulting in 



2c 



YJia 
1287r 



+ 



(A29) 



The " derivative" term 1^ can be easily calculated. The integration in the variables p and 
q again factorizes, and we have 



J// _ ^« , 
^ levr 

Putting all together, we finally arrive at the result 

Ci 4C2 



-ma + 



IGtt 



IGtt 



17 

21n2 



27rV2r2(i) 7ri/2r2(i; 



IGvr 



(A30) 



(Ina - 2.694) (A31) 



which is substituted in Eq. (p^). 



APPENDIX B: 

In this appendix we shall present the calculation of the integrals appearing in the elec- 
tromagnetic radiative corrections (Eq. (^)). 



{2txY (27r) 



g 2 _|_ ^2 



h+h 



(Bl) 



Integrating over the relative energy variables, the first term is rewritten in the form 
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e'^{M*f f d^p 1 f d?q 1 A-Kam^^)^ 



(27r)3 2 + ^2) 7 (27r)3 (w;(g))V2(p _ g)2 (^2 + ^2)2 



4my' J (27r)3«;(p)(p2 + ^2)2 
Using the same trick as in Eqs. (|Al| ) - (|A3|) , we can write 

' m'J'V 2n^Jo (l+p2)i/2(i + (i+p2)i/2)+---; ^1/2 1,^ ^)^--- ^^^^ 

In the calculations of I2 the term containing 2po'?o vanishes since it is odd in po and go- 
Thus one can write I2 = I2 + I'i where J2 and contain p^ and gg, respectively. 
I2 is ultraviolet divergent. Introducing dimensional regularization, we can write 



dT'p 



_(^2)2-n/2 



(277)' 



f + pf - ml^ (^(f - pf - m2^ - g )2 



167r2 V |g 



«;(g) - 




w{q) + 


1^1 



and 



^2 = ^(iV. + 2)^ (B5) 

/2 does not contain the ultraviolet divergence. Integrating over relative energy variables, 
we obtain 

jn _ /■ d?p d?q 1 1 47ram^0o 

'"47 (27r)3(27r)3u<p)(«;(g))i/2(p'_g')2(p'2 + ^2)(^2 + y) I ^) 

It is easy to see that I2 leads to the modification of the vr+vr" decay width in the order 
ct^lna and thus can be safely neglected. The final result reads as (cf. with Eq. (^51) ) 
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Table I. Corrections to the tt+tt atom decay width 



Effect 


Value 


Correction (in %) 


Strong 


5s 


-0.22 


Coulomb photon exchange 


Sc 


+1.55 


Mass shift 


5m 


+2.99 ±0.77 


Electromagnetic radiative 


5em 


+1.73 + 2.31 


Vacuum polarization 


5vac 


+0.19 


Finite size 


5f 


-0.16 


Total 


5tot 


+6.1 + 3.1 
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FIGURE CAPTIONS 



Fig. 1. Diagrammatic representation of the Bethe-Salpeter equation for the tt+tt" atom 
w.f. Initial equation (a) through the redefinition of the kernel (b) and the w.f. is transformed 
into the equation (c). The new kernel V contains the self-energy insertions in the outgoing 
external lines only (see (b)). 

Fig. 2. Matrix element describing the "residual" photon exchange. The imaginary part 
of the diagram (a) vanishes at the bound-state energy. Diagram (d) which is of the second 
order in the "strong" amplitude T12, is neglected in the present approximation. Dashed line 
denotes the virtual photon propagator, and dots correspond to the instantaneous virtual 
photon exchange. 

Fig. 3. Matrix element corresponding to the self-energy correction in the external pion legs 
(Eq. (H)). 

Fig. 4. Vacuum polarization correction. 
Fig. 5. Vertex correction. 

Fig. 6. Cancellation of the divergences which are present in the expression for the vertex 
correction. The divergence in the vertex subdiagram is cancelled by the counterterm depicted 
in (b) whereas the remaining overall divergence in diagram (a) is cancelled by the counterter 
given in (c). 
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